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ABSTRACT 



A computer ./var game is developed to measure the effect of false 
contacts on the probability of detecting a submarine. The variables 
are the probability of correctly classifying a non-submarine contact, 
the probability of correctly classifying a submarine contact, and the 
false contact density. A scenario is developed to focus on the false 
contact problem while holding other ASW variables constant. It is 
concluded from the output of the game that the effect of false contacts 
is deeply embedded in the interrelationships between units. 
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I. INTRODUCTION 



Considerable study effort has been expended on the problem of 
determining the probability of detecting a submarine in various Anti- 
Submarine Warfare (ASW) environments and tactical situations. In the 
phase of search and detection, classification of sonar contacts has been 
recognized as a serious complication. This complication develops 
because there are many substances in the sea that reflect a sonar 
signal in much the same way as a submarine. Thus, a sonar contact 
that is in fact a submarine may be classified as non-submarine while 
a non-submarine or false contact may be classified as a submarine. 

In the first case the searching units continue looking for a submarine 
that has already been found while in the latter case the time available 
for the search is expended prosecuting a false contact. 

Much has been done to improve the classification capability of 
ASW forces but the problem continues to be a serious one. In many 
ASW studies however, classification of sonar contacts is not considered. 
In Danskin's study [1] non-submarine contacts are not introduced into 
the problem and consequently it is neither necessary to revisit an area 
once it has been searched nor is time wasted prosecuting false contacts. 
Hammon [2] takes a typical war-gaming approach to the problem by 
assuming that for each helicopter a false contact is generated on any 
dip with a specified probability, and if a false contact is developed it 
is prosecuted by that one helicopter for a random period of time before 
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it is correctly classified. It is implicitly assumed that this in no way 
influences the search for the submarine by the other helicopters. 

The conclusions of both studies are valid if there exists a suitable 
degradation factor in the probability of detection that accounts for the 
effect of false contacts on the operational relationship among searching 
units. Unfortunately, the probabilities of correct classification are 
neither known with any degree of certainty nor are they static. In the 
Major Fleet Escort Study [3] it was assumed that a learning curve 
existed and that the probability of making an attack on a false contact, 
as a consequence of an incorrect classification, was a decreasing 
function of time. 

The purpose of this thesis is to examine the effect of false 
contacts on the probability of detecting a submarine in order to develop 
the relationship between the interactions among units and the probabil- 
ities of correctly classifying sonar contacts. This thesis will develop 
probability of detection curves as a function of the probabilities of 
correct classification for all values between zero and one. 
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II. SCENARIO AND ASSUMPTIONS 



In order to isolate the effect of false contacts a scenario was 
developed which held constant everything except those variables which 
bear on the classification problem. Although this creates a rather 
fictitious state of nature it does focus on the problem in its simpliest 
form much the same as a partial derivative represents the rate of change 
of a complicated function with respect to only one variable. 

The scenario is initiated by a submarine being detected by an 
outside agent. Neither the type of initial detection nor the length of 
time that the submarine is tracked is important. The only assumption 
necessary is that the contacting agent is certain that the contact is a 
submarine and can determine the position of the contact accurately. 

The contacting agent calls for helicopter assistance and then loses 
contact establishing the location of the last known position of the 
submarine (DATUM). The submarine assumes helicopters have been 
called in; dives to its best depth to avoid helicopter sonar detection, 
and proceeds to clear DATUM at a constant, randomly selected, speed 
less than or equal to its maximum speed available. The helicopters 
arrive at DATUM at some time late and commence the search. They have 
a specified time on station available for the search and if the submarine 
is not found by that time the search is discontinued. It is assumed that 
the helicopters know the maximum speed available to the submarine, 
but they cannot deduce either speed or course intelligence from the 
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tactical situation. They therefore commence a random dip search within 
the submarine's farthest -on circle allowing only a specified fraction of 
their detection range capability to overlap with any other unit's search 
area or the farthest-on circle. For its part the submarine can gain no 
intelligence from the active sonar transmissions of the helicopters and 
thus elects to take no evasive action. It is also assumed that the 
original contacting agent offers no further assistance to the helicopters. 

The helicopters have a "cookie cutter sonar" which always detects 
a contact that is within range and never detects a contact that is not 
within range. When a helicopter detects a contact the classification 
process begins. 

It is assumed that during the classification process there exists 
a constant probability that a valid contact will be correctly classified. 
There also exists a probability of correctly classifying a non-submarine 
contact however, in this case, it is assumed that the decision process 
in classifying a contact is influenced by previous classifications of the 
same contact. Thus it is assumed that a contact that has already been 
incorrectly classified is more likely to be incorrectly classified again 
than is a contact that is being classified for the first time. Alternatively 
the probability of generating a false contact prosecution is less than the 
probability of incorrectly classifying a contact given that it has already 
been incorrectly classified. 

If the submarine is detected and correctly classified the search is 
successful. If the submarine is incorrectly classified the helicopter 
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will continue searching for other contacts . If no others are found and if 



the helicopter has not reached the end of its dip cycle then the submarine 
contact will be re-evaluated. 

If a false contact is detected and correctly classified it is marked, 
so that no other unit will subsequently consider it, and the search is 
continued for the submarine. If the false contact is incorrectly 
classified the helicopter goes into the track mode and immediately re- 
evaluates the contact. The effect of two successive classification 
opportunities is to raise the probability of correctly classifying the 
initial contact in line with the previous discussion. If the contact is 
again incorrectly classified then a false contact prosecution is generated. 

If the number of helicopters prosecuting a false contact is less 
than a specified number and if another helicopter is available then it 
is called in to assist. When it arrives all units on the scene re-evaluate 
the contact. If the joint classification is correct the contact is scrubbed 
and all units return to the submarine search. Otherwise, all units 
continue to track the contact. Helicopters never lose contact on the 
false contact. This process is repeated incrementing the number of 
helicopters on the scene until either the contact is correctly classified 
or time on station is exceeded. The overall probability of detection is 
then computed as the number of times the search is successful divided 
by the number of times the scenario' is run. 
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HI. DESCRIPTION OF THE MODEL 



A Monte Carlo event store computer war game, named "Helicopter 
DATUM Search" (HSDS) , was developed to implement the scenario and 
assumptions. The principals of the game are the helicopters, submarine, 
and the false contacts. The helicopters have attributes of location and 
status while the submarine and and false contacts have attributes of 
location, velocity and status. HSDS can handle variations in tactics 
of either the submarine or helicopters. Up to 6 helicopters and up to 
1000 false contacts can be introduced with any combination of time 
late, time on station, and probabilities of correct classification. This 
provides the capability within HSDS of varying the density of false 
contacts and the size of the initial area that the submarine can be in 
when the helicopters first arrive at DATUM. 

The logic of the game is written in FORTRAN IV for the IBM 360 
computer and several support routines are written in OS/3 60 Assembler 
Language. The entire game is modularized with each function programmed 
as a separate subroutine to permit ease in modification for future study. 
The entire program requires less than 100K bytes of core storage and one 
run of 1000 replications with 4 helicopters and 250 false contacts requires 
slightly less than 5 minutes computing time. 
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A. EXECUTIVE ROUTINE 

This section contains a description of the executive routine, or 
MAIN program, followed by the subroutines that support the MAIN program. 
These routines are written in FORTRAN IV and listed in Appendices B and 
D. A diagram of the flow between the MAIN program and the Events is 
contained in Figure 1 . 

The MAIN program first reads the input data. The salient parts 
of the input data are then printed and the run and first game are 
initialized. The helicopters start at DATUM. The submarine's speed is 
determined by a random number multiplied by the submarine's maximum 
speed. False contacts are distributed uniformly within an area bounded 
by the submarine's farthest-on circle at the end of time on station. A 
control loop is then entered taking the top event off the calendar and 
checking to see if time on station has been exceeded. If so , a branch 
is made to Subroutine OVER. Otherwise, the appropriate routine is 
called and upon normal return from that routine the program loops back 
to get another event from the calendar. Upon normal return from sub- 
routine OVER the program loops back to initialize the next game. If 
this was the last replication of the run the results of the current run 
are printed and the MAIN program loops back to read another input card 
for the next run. If all cards have been read, and thus all runs are 
complete, the program terminates. 

Subroutine OVER computes statistics for each game. The total 
number of replications of the game is divided into two groups. At the 
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GENERAL PROGRAM FLOW CHART 



end of the first group of replications the random number generator is reset 



and the sign of the argument for the random number generator is changed. 
This causes the next stream of random numbers to be the compliment of 
the previous stream. The resulting statistics are assumed to have the 
Antithetic properties discussed in Hillier and Lieberman [4] . This is a 
technique which uses negative correlation between groups to produce a 
negative covariance thereby reducing the variance of the average of the 
groups. A special return to the MAIN program is made if the last 
replication of all groups has been played to indicate the end of the run. 

Subroutine SETUP computes the time interval for each branch to 
the SEE event and schedules the SEE event on the calendar. 

B. EVENTS 

This section contains a description of the event routines followed 
by the subroutines that support the event routines. Each of these 
routines are written in FORTRAN IV and listed in Appendices B and D. 

Event BRKDIP calculates the next dip position of a specified 
helicopter by using two random numbers to compute the X and Y coordinates 
of the new DIP station. If this position is within the farthest-on circle 
of the submarine and does not overlap another helicopter's search area 
by more than the specified amount, then a DIP event is scheduled for 
that position after calculating the time to station. Helicopters transit 

v 

between DIP stations at 90 knots and require 3 minutes for lowering 
and raising the sonar dome. If either of the above tests fail a new 
position is calculated. Note that an infinite loop can be generated here 
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if the sum of the helicopter's search areas is greater than the area 
within the submarine's farthest-on circle. 

Event DIP handles the helicopter in the search phase. It first 
schedules a BRKDIP event and then branches to subroutine DET for 
detections. A normal return to the MAIN program is made if there were 
no detections. If a detection was made a branch is made to subroutine 
CLASS to determine the classification. If the submarine was the contact 
detected and it is correctly classified a special return is made to the 
MAIN program to indicate the search has been completed successfully. 

If the submarine was incorrectly classified a RELOOK event is scheduled. 
RELOOK is an entry point in the DIP routine just before the call to sub- 
routine CLASS. If a non-submarine contact was detected and correctly 
classified the non-submarine's status is changed to inactive and a 
branch is made to subroutine DET to determine if any more detections 
are possible on this dip. If a non-submarine was detected and incorrectly 
classified a TRACK event is scheduled for the current game time. 

Subroutine CLASS determines the classification of a contact, 
submarine or non-submarine, by comparing the input probability of 
correct classification with a random number. 

Subroutine DET computes the index of the first contact that is 
within sonar range of the specified helicopter. Only active contacts 
are considered for detection. This routine gives a slight advantage to 
the helicopters since the submarine is always the first contact considered. 
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Subroutine UPDATE computes the position of all active contacts 



for the present game time. 

Event SEE prints the position and status of all helicopters and 
contacts at the present game time. 

Event TRACK handles the prosecution of false contacts. A branch 
is made to subroutine JOINT to determine the current classification. 

If this classification is correct, a branch is made to subroutine RESUME 
to discontinue prosecution of the contact. If the current classification 
is incorrect another TRACK event is scheduled and the number of 
helicopters on the scene is determined. If the number on the scene is 
less than the maximum number allowed and if there is another helicopter 
enroute to a new DIP station or in a DIP but not prosecuting a contact, 
that helicopter is sent to the scene. The position of all helicopters on 
the scene is then updated to the false contact's present position. 

Subroutine JOINT determines the composite classification of the 
contact by all of the units on the scene prosecuting the false contact. 
The criterion for the classification is an average value computed as 
follows: 





NOS(K) 


SUM = 1/NOS (K) 


x x(j) 

j=i 


where X (j) 


= 1 for RN(IX) < PCNS 


X(j) 


= 0 for RN (COUPONS 


NOS(K) 


= number of helicopters prosecuting contact K 


PCNS 


= probability of correctly classifying a false contact 


RN (DC) 


= random number from stream IX 
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The classification is correct if SUM is greater than PCNS and incorrect 
otherwise . 

Subroutine RESUME calls BRKDIP for each helicopter that has been 
prosecuting a false contact and then sets the inactive code for that 
contact. 

C . SUPPORT ROUTINES 

The following routines, written in OS/3 60 Assembler Language, 
are general to any simulation and are not included in Appendices B and D. 

RN is an additive random number generator which can produce a 
minimum of 16 unique random number streams. If the argument to the 
function is zero the generator is initialized and no number is returned. 

A positive argument produces a random number from the specified stream 
of numbers while a negative argument produces the compliment random 
number from the same stream. Two generators are used in HSDS; one 
for movement and position variables, and the other to determine 
classifications . 

Subroutines INT, TNE, SNE, and REMOVE comprise the event 
calendar package. INT initializes the calendar and must be called 
before the calendar is used. TNE takes the top event off the calendar. 
SNE stores an event on the calendar in its proper time sequence. 

REMOVE nullifies all events on the calendar which have a specified 
value for a specified attribute. In HSDS events are removed from the 
calendar as they are invalidated rather than checking each event as it 
is returned by TNE. 
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Subroutine CALFUL, called when an error condition develops, 



terminates the program and produces a core dump. 

Subroutine MOMENT calculates the mean and standard deviation 
of an input group of data . 



17 



IV. VALIDATION AND ANALYSIS 



The first objective of this thesis was to measure the effect of the 
probability of correctly classifying a non-submarine on the overall 
probability of detection. Within the general scenario of Chapter II, a 
set of conditions, designed the base conditions, was chosen as follows: 



Number of helicopters 4 

Time late to DATUM 15 min. 

Time on station 60 min. 

Expected sonar range 3000 yards 

Maximum submarine speed 20 kts . 

Number of false contacts 250 



Probability of correctly classifying a submarine contact 0.8 
This corresponds to a false contact density of 0 . 1 contacts per square 
mile; a probability of detecting a false contact of 0.7 for any helicopter 
on any dip, and a probability of detecting the submarine on the first dip 
after reaching DATUM of 0.053 for any helicopter. This last probability 
is computed using time late plus the expected value of 4 . 5 minutes for 
the time from DATUM until the helicopter is in the first dip and searching 
for the submarine to determine the submarine's farthest-on circle. The 
probability of detection is then the area of the helicopter's sonar search 
divided by the area within the submarine's farthest-on circle. These 
conditions were chosen because other studies have considered the same 
or similar conditions. 
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A. RESULTS AND VALIDATION 

Before a production run was made for each parameter, a pilot study 
was conducted to determine the number of runs that would produce 
statistics with an acceptable variance in the minimum amount of computing 
time. The number of runs described in the following sections are a result 
of such a pilot study. 

The probability of correctly classifying a non-submarine contact 
(PCNS) was incremented from 1.0 to 0.0 and at each point 2 000 repli- 
cations of the game were run and the probability of detection calculated. 
The manner in which helicopters interacted was considered to be the 
maximum number of helicopters that would prosecute a contact (MAXNOS) . 
As a result of these runs the probability of detection as a function of 
PCNS and MAXNOS is tabulated in Table 1 . Additionally since a false 
contact must be incorrectly classified twice in succession to generate 
a false contact prosecution, PNFC is defined as the probability of NOT 
generating a false contact prosecution. PNFC is computed as: 

PNFC = 1.0 - (1.0 - PCNS) 2 
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MAXNOS 



PCNS 


PNFC 


1 


2 


3 


4 


1.0 


1.0 


.371 


.371 


.371 


.371 


0.8 


0.96 


.394 


.373 


.372 


.359 


0.6 


0.84 


.390 


.346 


.348 


.350 


0.4 


0.64 


.371 


.307 


.315 


.322 


0.2 


0.36 


.342 


.266 


.259 


.249 


0.1 


0.19 


.316 


.221 


.211 


.218 


0.05 


0.0975 


.303 


.196 


.179 


.173 


0.0 


0.0 


.256 


.162 


.160 


.149 



The Probability of Detection 
TABLE I 

For each value of the maximum number of helicopters that will prosecute 
a submarine, the probability of detection is well approximated by a 
linear function in PNFC . A representative linear fit of the data is 
plotted in Figure 2. The linearity assumption between the probability 
of detection and PNFC appears reasonable since, with independent 
random dipping, the probability of detection should be degraded only by 
time lost in prosecuting false contacts. This is in contrast with a search 
plan where helicopters work together and the effectiveness of the search 
plan is dependent on the completion of a particular pattern. The probability 
of detection when PCNS is 1.0 is also reasonable. In LCDR Buck's 
study [5], which does not consider false contacts, the probability of 
detection for two helicopters against a closing submarine under the same 
conditions was computed to be 0.37. In the present analysis there are 
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0.38 



Pd 



l 



i 



i 





0.0 0.25 0.50 0.75 1.0 

PNFC 

PROBABILITY OF DETECTION (Pd) 



vs 

PROBABILITY OF NOT GENERATING 
A FALSE CONTACT (PNFC) 



FIGURE 2 
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four helicopters searching for a submarine that can depart DATUM in 
any direction and the probabilities of detection are in close agreement. 

B . ANALYSIS 

The results of Table I demonstrate how the effect of false contacts 
varies with the extent of helicopter interaction. Figure 3 is a plot of 
the least square fit of the data taken of the probability of detection as 
a function of PNFC for the maximum number of helicopters that will 
prosecute a contact (MAXNOS) . The gap at PNFC equals 1.0 is the 
result of the least square approximation. 

Because of the simplicity of the scenario and tactics, it is possible 
to compute the probability of detection for all of the helicopters on the 
first dip analytically. Consider PNFC equal zero so that any helicopter 
detecting a false contact is lost to the search for the remaining time 
left on station. For MAXNOS of 1 the probability of detection for all of 
the helicopters on the first dip is equal to the sum of the probabilities 
that each helicopter detects the submarine on its first dip multiplied 
by the probability that the submarine is correctly classified. This is 
so since the actions of each helicopter are independent. If a false 
contact is detected and hence incorrectly classified by a helicopter, 
it is prosecuted by that one helicopter while the other helicopters 
continue the search unaffected. For the base conditions the overall 
probability of detection is approximately 
4 x 0.053 x 0.8 = 0.17 
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When MAXNOS is greater than 1 the actions of each helicopter 
are no longer independent. The first helicopter to dip detects the 
submarine with probability 0.053 and detects a false contact and in- 
correctly classifies it as submarine with probability 0.7. If the first 
helicopter detects and incorrectly classifies a false contact the next 
helicopter comes to the first helicopter's assistance rather than con- 
tinuing to its first dip station and thus never gets into the search for 
the submarine. The actions of each helicopter in turn depend on what 
has occurred to the previous helicopters. The overall probability of 
detecting the submarine for all of the helicopters is then the uncon- 
ditioned probability of the sum of each conditional probability which is 
calculated below. In these calculations it is assumed that on the first 
dip the time required for a joint classification exceeds the time to 
station for subsequent helicopters and therefore no more than one 
helicopter can come to the assistance of another helicopter. 



Let: 



A 



The event the submarine is correctly classified 
given that the submarine is detected 



B/j 



The event the submarine is detected by helicopter 
j given that helicopter j is not diverted 



C/1 



The event that helicopter 1 is not diverted which 
always happens 



C/2 



The event that helicopter 2 is not diverted which 
happens only if helicopter 1 does not detect and 
incorrectly classifies a false contact 
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C/3 = The event that helicopter 3 is not diverted which happens 

only if: 

1) helicopter 1 detects and incorrectly classifies a 
false contact thus diverting helicopter 2, or 

2) neither helicopters 1 nor 2 detect a false contact 

C/4 = The event helicopter 4 is not diverted which happens 

only if: 

1) helicopter 1 detects and incorrectly classifies a 
false contact diverting helicopter 2 and helicopter 
3 does not detect a false contact, or 

2) helicopter 1 does not detect a false contact and 
helicopter 2 does detect and incorrectly classifies 
a false contact diverting helicopter 3, or 

3) neither helicopters 1,2, nor 3 detect a false contact 



Then: 



Probability of A 


= 


PCS 




— 


0.8 


Probability of B/j 


= 


Probability of detection multiplied by 
the probability the submarine is not 
detected by a previous helicopter 




= 


0.053 x 0 . 954^ _1 


Probability of C/1 


= 


1.0 


C/2 


= 


0.3 


C/3 


= 


0.7 + 0.3 2 


C/4 


= 


0.7 x 0.3 + 0.3 x 0.7 + 0.3 3 


The overall probability 


of detection is: 



4 

n P(A) X P(B/j) x P(C/j) 
j=l 

= 0.101 

This is in comparison with the 0.17 probability of detection when there 
is no interaction between helicopters. 
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When PNFC is greater than zero the analysis is the same except 
that the conditional probabilities are further conditioned on the generation 
of a false contact prosecution given that a false contact has been detected 
by a previous helicopter. The difference in the probability of detection 
between MAXNOS of 1 and values greater than 1 approaches zero as PNFC 
approaches one. 

The total probability of detection is of course greater than the 
probability of detection on the first dip but an analytic computation for 
subsequent dips becomes quite involved since the single dip probability 
of detection decreases as the square of the submarine's maximum speed. 
What the analytic computation for the first dip shows is that the dis- 
proportionate degradiations in the probability of detection for MAXNOS 
of 1 , 2, and 4 is neither an accident nor a mistake. Thus even when 
the probability of detection is a simple linear function of the probability 
of prosecuting a false contact for any individual unit, the probability 
of detection for the group varies in a much more complicated manner. 

C. EFFECT OF OTHER VARIABLES 

Since the probability of detection is conditioned on the probability 
of correctly classifying a submarine contact (PCS) and the probability 
of detecting a false contact, or false contact density, it is necessary 
to determine how the probability of detection varies as a function of 
PCS and the false contact density. Although an argument can be made 
to the contrary, it is assumed here that these two phenomena are 
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independent and can be treated separately. No assumption is made that 



either PCNS and PCS or PCNS and the false contact density are independen 
PCNS was fixed at 0.4 and the probability of correctly classifying 
a submarine was incremented from 1.0 to 0.0 with 1000 replications 
made at each data point. The results are well approximated by a fourth 
degree least square polynomial as plotted in Figure 4 . The significance 
of this curve is both in its effect on the conditional probability of 
detection as previously discussed and on the design trade-offs in 
hardware. In many cases an increase in the range of a detector, and 
thus an increase in the conditional probability of detection of the sub- 
marine, is accompanied by a degradation in the classification clues. 

Such is the case when going to lower frequency sonars . Under 
conditions where the probability of correct classification is high, say 
greater than 0.6, a greater degradiation can be accepted in classification 
clues, or classification capability, than when the probability of correct 
classification is low, say less than 0.4. This is because the slope of 
the probability of detection curve is decreasing with increasing probability 
of correct classification. 

Figure 5 is a plot of the least square linear fit of the probability 
of detection as the false contact density was incremented from 0 . 0 to 
0.2 false contacts per square mile. PCNS and PCS were set at 0.4 and 
0.8 respectively and the maximum number of helicopters to prosecute 
a false contact was set at 4. Each data point is the result of 1000 
replications of the game. Least square polynomials of higher degree 
were investigated but provided little additional information. 
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PROBABILITY OF DETECTION (Pd) 
vs 

FALSE CONTACT DENSITY PER SQUARE MILE ( FCD) 



v. CONCLUSIONS AND RECOMMENDATIONS 



The analysis of Chapter IV indicates that the effect of false contacts 
is not obvious even under the simpliest of conditions. It is not un- 
reasonable to assume that as the tactics become more involved the effect 
of false contacts becomes more deeply embedded in the interrelationships 
between units. This is particularly true when different types of units 
with different sensors and classification capabilities operate together. 

The analytic formulation of such a problem is formidable is not impossible. 
However, if the interrelationships can be well defined such a problem 
can be formulated as a game. 

It is important to recognize that false contacts can play as large 
a part in the probability of detection of a submarine as any other 
parameter. It is therefore recommended that classification and the 
effect of false contacts be considered in future ASW studies. This could 
be accomplished by running a pilot study of the tactics under consideration 
in a simulation, such as the game in this thesis, in order to develop the 
appropriate degradation curves for those tactics and forces. It is further 
recommended that future ASW war games consider the interrelationships 
between units when developing the logic for false contact prosecutions. 
This is of particular importance where the game will be used to help in 
the evaluation of new equipments where trade-offs have been made in 
the classification capability. 
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There are several areas within this thesis which could be fruitful 



topics for further study. Current helicopter search plans could be 
implemented, by changing the BRKDIP routine, in order to study their 
sensitivity to classification errors. Subroutine JOINT could be modified 
to provide for the inclusion of different types of sensors to aid in the 
classification process when more than one unit is on the scene. Finally, 
a submarine maneuver event could be added to study the effect of 
submarine evasion. 
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APPENDIX A 



LIST OF VARIABLES 

A. INPUT VARIABLES 

The following is a list of input variables in the order that they are 
read into the program: 



IN 


Input data set number 


OUT 


Output data set number 


IRUN 


Number of this run 


NGRPS 


Number of groups for each run 


NREPS 


Number of replications for each group 


IX 


Argument for the random number generator 


NN 


Number of intervals in the first game that the 
positions of all units will be printed 


NHS 


Number of helicopters 


TL 


Time late 


DELT 


Duration of each dip 


TOS 


Time on station 


RS 


Maximum sonar range 


OVRLAP 


Amount of overlap permitted between two HS sonar 
search areas 


NUM 


Number of false contacts (The program increments 
this number by 1 to include the SS in all calculations) 


VMAX 


Maximum speed of the SS 


PCS 


Probability of correctly classifying a submarine 
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PCNS 



Probability of correctly classifying a non-submarine 



MAXNOS Maximum number of HS that will prosecute a false 
contact 

B. INTERNAL VARIABLES 

The following is a list of variables internal to the program as they 
appear in common: 



XS(J) 


X position of contact J 


YS(J) 


Y position of contact J 


VXS(J) 


X velocity component of contact J 


VYS(J) 


Y velocity component of contact J 


ISTAT(J) 


Status of contact J 
0 => Active, 1 => Inactive 


TS (J) 


Time contact J was last updated 


NOS(J) 


Number of HS prosecuting contact J 


MAX 


Minimum of (NHS , MAXNOS) 


XH(J) 


X position of helicopter J 


YH(J) 


Y position of helicopter J 


ICODE(J) 


Status of helicopter J 

-1 => Airborne enroute to next dip 

0 => In dip searching 

N > 0 => prosecuting contact N 


TP 


Current game time 


IDUM 


Dummy variable for arguments to SNE 


NS 


Number of times the submarine was detected 
in this group 


TAVG(N) 


Time to detection for detection N 


NGAME 


Counter for the number of games played 
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C . OUTPUT VARIABLES 



The following is a list of variables output by the program in the 
order that they are printed: 



XBAR 


Average probability of detection 


SD 


Standard deviation of the probability of detection 


AVG(J) 


Probability of detection for group J 


TBAR 


Average time to detection given that a detection 
was made 


TSD 


Standard deviation for time to detection 


NC 


Total number of times the submarine was detected 
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APPENDIX B 



FLOW CHARTS 



The following is a list of symbols and their meanings used in the 
flow charts of HSDS. 






Entry or exit point 



Input/Output operation 



Decision point 



Calculation or Instruction 



Set of calculations expanded 
on a following flow chart 



/ 



\ 



( 



Store an event on the calendar 



Multiple branch 
Computer GOTO 




Off-page connector 
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APPENDIX C 



INPUT AND OUTPUT 



A. 


INPUT FORMAT 










CARD 

NO. 


VARIABLE 

NAME 


CARD 

COLUMN 


FORMAT 


RANGE 


UNITS 


1 


IN 


] 


11 






1 


OUT 


2 


11 






2 


IRUN 


1-5 


15 






2 


NGRPS 


6-10 


15 


NGRPS 2 




2 


NREPS 


11-15 


15 


NREPS <.99999 




2 


IX 


16-20 


15 


1 < IX < 16 




2 


NN 


21-25 


15 






3 


NHS 


1-10 


110 


1 < NHS £6 




3 


TL 


11-20 


F 




minutes 


3 


DELT 


21-30 


F 




minutes 


3 


TOS 


31-40 


F 




minutes 


3 


RS 


41-50 


F 




yards 


3 


OVRLAP 


51-60 


F 


1.0£OVRLAP£2.0 
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CARD 

NO. 


VARIABLE 

NAME 


CARD 

COLUMN 


FORMAT 


RANGE 


UNITS 


4 


NUM 


1-10 


no 


0<NUM1999 




4 


VMAX 


1-20 


F 




knots 


4 


PCS 


21-30 


F 


O.OiPCS^l.O 




4 


PC NS 


31-40 


F 


O.O^PCNS <1.0 




4 


MAXNOS 


41-45 


15 


MAXNOS >0 




NOTE: 


F denotes a 


real variable 


requiring 10 


card columns and a 





decimal point. 

Ixx denotes an integer variable right-justified in xx card 
columns . 
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B 



SAMPLE OUTPUT 
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noon oooo oooo oooo oono oooo 



APPENDIX 0 



PROGRAM LISTING 



CCNNCN /SS/ XS( 1000) ,YS(1000) t VXS( 1000 ) ,VYS( 1000) , 

* I ST A T ( 1000 ) , TS( 1000) , VMA X ,N J St 1 COO ) ,NUM,MAX 

CCMMCN /MS/ XH( 6 ) , YH( C ) , IC0DE ( 6 ) ,RS ,CELT , NFS ,TL ,TGS 

CCMMCN /UT IL/T P , IX, PCS, PCNS , I CUR, NS , AVG( 2 ) , T AVG( 1000) , 
*NC,TBAR, TSC,NGAME,NREPS ,NGRPS, I CUT,CVRLAP,XBAR,SD 

READ INPUT CATA 



PEAD(A, 199 ) IN, ICUT 

READ( IN, IOC) IRUN,NGRPS,NREPS,IX » NN 
RE AD IIN, 10 1 )NHS,TL,DELT,TOS,RS,CVRLAP 
R S =P S/2000.0 

7CC READIIN, 102, END =90C)NUR,VRAX, PCS, PCNS, MAXNCS 
MX=MINO(MAXNUS,NHS ) 

WRITE INPUT PARAMETERS FUR THIS RUN 



WF I T E ( IGUT.200) IRUN 
I F LN= I PUN + 1 
PPS=RS*20CC.C 

WRITE ( ICUT ,201 ) NR EPS, I X,NHS,TL ,TCS,RRS 
WRITE! I CUT, 202 ) VMAX ,NUM , PCS ,PCNS , MAX 

IN IT IALIZE this run 



NMAX=NP EPS*NGRPS 
NCAME=1 
NLM=NUM+ 1 
IE UM=- 9999999 
NC = 0 

S LM=RN 1 ( C) 

SLM = RN2 ( 0 ) 

N S = 0 

CALL INT 
CALL SETUP ( NN ) 

CCTC 501 

INITIALIZE THIS GAME 



5C0 CALL INT 
501 T P = T L 

SET MS STARTING DIP STATIONS 



CC 301 1 = 1, NhS 
XM I ) =0 .0 
Yh(I )=0.c 
3 C 1 CCNTINUE 

CC 3C2 1=1 , NHS 
CALL BRKDI P( I ) 

302 CCNTINUE 

SFT SS CCUPSE SPEEC E STATUS 



VX$< 1) =PN1 ( I X) *VMAX 
VY S ( 1 ) =0 .0 

XS ( 1 )=VXS( 1 > * T L/60 . 0 
YS ( 1 ) = 0.C 
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moo moooooo oooo nnnnmn noon noon 



1 S ( 1 ) = T L 
I 5 T A T ( i |=0 

INITIALIZE HLS[ CONTACTS 



IF (MiM-1) 401,401,400 

'.CO P = TCS* VPAX/60. 0 
EC 303 1=2, NUM 
X S ( I ) = «M< IX )*2.0*R-R 
V S < I )=KN 1 ( I X) *2.0*R-R 
VXS < I ) =0.0 
VYS ! I ) =2.0 
NC S ! I ) = C 
ISTAT ( I ) =0 
T S ( I ) = T l 
3 C ? CONTINUE 

TAKt NEXT tVENT ANC CF'ECK TI^p 



<*C1 CALL TNf!TTP t ISUH» I UNIT » [ V AL » 6 y 50 ) 

IT (TTP-TP) 701,702,702 
701 CALL CALfUL ( 15 I 
7 C 2 T F = T TP 

IF (TP+TL-TGS) 402,^02,403 

GC TL PROPER EVENT ROUTINE, EXECUTE EVFNT 





TFFN 


return for next event 


A C 2 


.GC TC 


( 6 Cl , 602, 6C 3, 6 04 ,60 5 ,606 ,60 7,6 08 ) , I SUP 


tCl 


CALL 


bRKCIP! IUN I T ) 




CC TO 


401 


6C2 


CALL 


DIP! I UN I T ,640 ? ) 




GC TO 


40 I 


c03 


CALL 


GET ( IUN IT , I V A L ) 




GC TC 


401 


6 04 


CALL 


TRACK ( IVAL ) 




GLTC 


401 


r C 5 


CALL 


RE LOCK! I UN I T , I VAL ,6403 ) 




C C TC 


401 


6 C 6 


CALL 


SEE 




GL TO 


401 


6 C 7 


CLTC 


40 1 


6 C 3 


CC T C 


401 




THIS 


game IS COMPLETE 



4C2 CALL CVER! £404 ) 
NGAME=NGAME+1 
GC TC 500 

T I I S RUN I S COMPLETE 



WRITE RESULTS AND CC BACK F C R ANCTHER RUNi 



4C4 WRITE! I CUT, 203) X BAR, SC 

*RITE! I OUT ,204M AVO( I ) , 1=1 ,NGRPS) 
WRITE! I CUT ».?05)TBAR,TSC,NC 
GC TO 700 

ALL RUNS COMPLETE 



5CC S T C P 

SFF CALL CALFUL! 10 ) 
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noon rrooo rrooo 



199 FCRMATI2I1) 

ICO FCRMAT (515 ) 

101 FORMAT! 1 10 i 5F10.0) 

2 CO FORMAT (/// 1H1 i^Ox’, 30HHS DATUM SEARCH *** RUN NUMBER, 
*15) 

201 FCPMAT(lH0,20X,34HTFE INPUT PARAMETERS PCP THIS RUN , 

^CHTHE’PANrCM 8 NUMBER SECUENCe’iS , I 2 , // 1 OX , 9HTHE R t ARE 
*»I3t32H HELICOPTERS WITH A TIME LATE OF,F6.2,PH MINUTE 

4 4 F AND * A TIME CN STATICN GF,F6.l,8H M I NUTE S , / /I OX , 
*22FEXPECTFD SONAR RANGE =,FS.2,6H YARDS,/! 

202 FORMAT ( /10X , 30HTHE MAXIMUM SUBMAR I NE SP EEC l S , F5 . 1 » 

* 6 F KNOT S ,/ / 10 X , 10 HTHERF AKn ,I3,15H FALSE CuNTACTS,// 

* 10X , 42HTHE PRUBAdlLITY CF CORRECTLY CLASSIFING A , 

* 20HSUBMAR I NE C CN T AC T I S , F 5 . 2 , / / 10X , 1AHTH E FRCBAEI LITY 
*50FOF CORRECTLY CLASSIFING A NON- SUBMAR I NE CONTACT IS, 
*F5.2»//10X»7HAT MOST,I2,27H HELICOPTERS WILL PROSECUTE 

* 1 OF A CONTACT) 

203 FCPMiAT I///20X, 19HRESULTS OF THIS RUN.//10X, 

* 22HTHE PROBABILITY CF DETECTIUN W AS , F 8. 5 , / / 1 OX , 

* 30FW 1 TH ONE STANDARD DEVIATION OF.F8.5/) 

204 FCRMAT ( /10X, 31HTHE PPG FAB I L IT I ES BY GROUP WERE,F8.5, 
*4F A ND , T B. 5/) 

-05 FCRMAT (/ 10X, 33HTHE AVERAGE TIME TO DETECTION WAS,F10.5 
*/ / 10 X * 30HW I T H ONE STANDARD DEVIATION OF , F 1 C , 5 , / / 1 OX , 

* 2 5FBASED UN A SAMPLE SIZE OF, 14) 

ENC 



SLBRCUTINE ORKDIP(N) 

CCMMCN /SS/ XS ( 1000 ) ,YS ( 1000 ) ,VXS( 1000 ) ,VYS( 1000 ) , 

♦ I ST ATI 1C00 ) , TSIIOOO) , VMAX ,NOS( 1000 ) ,NUM,MAX 
CCMMCN /HS/ XH(6),YF(6), ICODE ( 6 ) , RS , DEL T , NHS , TL , TOS 
CCMMCN /UT IL/TP , IX, PCS, PCNS, IDUM, NS, AVG( 2 ) , TAVGI 1000) , 
*NC ,TBAR ,TSD,NGAME,NREPS,NGRPS,ICUT,CVPLAP,X3AR,SD 

CANCEL ALL FUTURE EVENTS AND SET AIRBORNE COCE 



CALL P EMCV F 1 2 , N ) 

I C CDF ( N ) =- 1 

CCMPLTE ANC TEST NEXT DIP STATICN 



R=VMAX*TP/H0.0 
5CC X = RN 1 ( I X )*2. 0*R-R 
Y = R N 1 ( IX )*2.0*R-R 

IF ( ( X* X+Y*Y )- ( R-RS )** 2 ) 400,400,500 
4CC CC 3C0 1=1, NHS 

R 1= ( X-XF ( N ) )’**2+(Y-YH(N))**2 
IF (SGPT(Rl) - C VR L A P* R S ) 500,500,300 
300 CONTINUE 

COMPUTE TIME TC STATICN C STORE DIP EVENT 



TCE=TP+SCRT( (X-XH(N) )*v2+(Y-YH( N) ) **2 ) *0 . 667+3.0 
CALL SNE(TCE,2,N, ICUM, £999) 

XF ( N ) = X 
YF (N )=Y 
RETURN 

999 CALL CALFUL(l) 

END 
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S LOR CUT INF CLASS(JfK) 

C 

CCNWCN /SS/ XS < 1000 ), YS ( 1000 ), VXS( 1000 ), VYS( 1C0C ) , 

* I STAT( 1000 ) , TS( 1000) ,VWAX,N0S< 1000 I , N U M , WAX 
CLNF’CN /HS/ XH< 6 ) , Yh< 6 ) , ICOOEI 6 ) ,PS ,DELT ,NHS ,TL ,TCS 
CCWMCN /UT IL/T P , IX , PCS, PCNS, I HUP , MS , AVG ( 2 > , T AVG ( 1000) , 
*NC , TPAK , TSD.NGAWE , N RT P S , NGRP S , I CUT ,CVRLAP,XBAR, SC 

J IS THE SUBSCRIPT CF THE CONTACT 

K IS THE CLASSIFICATION CODE 

IF (J-l) AC2,ACC,AC2 
ACC K= 1 

IF ( RN2 ( I X ) - °C S ) 5 C C » 5 00 » AO l 
A C 1 K = 2 

C-C TO 500 
K = 2 

IF (RN2( IX J-PCNS) 5CC , 500 , AO T 
A 0 2 K = A 
5 C C RETURN 
ENC 



SLPKCUT INF CET(J,K) 

CCWNCN /SS/ XS ( 1000 ), YS ( 1000 ), VXS( 1000) ,VYS( 1CC0) , 

* I STATI 1 COC ) , TS < 1000 ) ,V M AX ,NOS ( 1000 ) , NUM, WAX 
CLWWLN /HS/ XH(F) ,YH(6) , ICC«Jn(t») , R S , I) t L T , NH S , T L , TGS 
CCNWCN /UTIL/TP,IX,PCS,PCNS, I DUN , NS , AVG ( ? ) ,TAVG( 1000) , 
»NC ,TBAR ,TSC tNGAPfc , NPEPSt NGRPSt I CUT , CVRLAP , X 3 A R , SD 
CALL UPDATE 

SET K = THE SUBSCRIPT CE THE FIRST CONTACT WITHIN 
RANGE UP = 0 IF NO CONTACT IS WITHIN P.ANGt 

DC 2 CO 1=1, NUM 
IF ( I S T AT ( I ) - 1 ) AOC, 30C,A0C 
ACC T ENP=<XEi<J)-XS<l))*«? + (YH(J)-YS<I)>**2 
IF ( TEM.P-RS-RS ) AO 1 , A01 , 300 
A 0 1 K = I 

GCTC 500 
2CC CONTINUE 
K — C 

5CC RETURN 
END 



SUBROUTINE D I P ( N , * ) 

CCWMCN /SS/ XS < 1000 ), YS ( 1000 ) ,VXS< 1000) , VYS( 1000) , 

* I STATI 1000 ) , T S ( 1000) , V WAX , NGS (1000) , NUM, WAX 
CCNMUN /HS/ XH(fc) ,Yh<6> , ICCDEI6) , K S , D£ L T , NH S , T L , TCS 
CCWWLN /UTIL/TP,IX,PCS,PCNS, I DUN , NS , AVG ( 2 ) ,TAVG( 1000) , 

* NC , TEAR ,TSD,NGAWE,NREPS,NGRPS, I CUT, CVRLAP , XBAR, SD 

SET DIP CUCE, STURE NEXT 3RKC IP EVENT 
i CHECK FCR DETECTIONS 

I CC C E ( N ) =0 
TC=TP+DELT 

CALL SME <TC , 1 , N , I DUN, 0559 ) 

5CC CALL DETIN, K) 

IF (K) A02 , AC2 , ACC 

ENTRY RFLPCKIN.K, *) 

CLASSIFY CCMACT 
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4CC C/'LL CLASS (K ,L ) 

GCTO(fcC0ffcCl»6C2*6O3) , L 

THE SS HAS BEEN UETECTEC & CORRECTLY CLASSIFIED 

6CC I S T A T ( 1) = I 
RETURN I 

THE FALSE CONTACT HAS 8EEN CORRECTLY CLASSIFIED. 
CELETE IT E CHECK FCR MORE CONTACTS 



601 I S T AT ( K ) = 1 
GLTC 600 

THE SS HAS <JEEN INCORRECTLY CLASSIFIED. 
SCHEDULE A RELCCK EVENT 



602 



T CE = TP + 1.0 

CALL SNEITCE, 5, N, 1,6999) 
GCTO A 0 2 



THE FALSE CCNT ACT HAS SEEN INCORRECTLY CLASSIFIED. 
SCHEDULE A TRACK EVENT FOR THE CURRENT GAME TIRE 



603 NC S I K ) = 1 

ICCDEI N ) =K 
CALL R EMOV E I 2 » N ) 

CALL SNE (T F , 4 » I CUM » K, 6999 ) 
A C 2 RETURN 
999 CALL CALFULI2) 

ENC 



SUBROUTINE JUINT(K.L) 

CCMMCN /SS/ XSI 100C) ,YS< 1000 ) ,VXS( 1000) ,VYS< 1000 ) , 

’I ST ATI 1000 I ,TS< 1000 ) , VM A X ,N0 S I 1 000 ) ,NUM,MAX 
CCMMCN /HS/ XH(6) » YF I 6 I » IC0DEI6 ),RS»DELT , NHS»TL»TOS 
COMMON /LTIL/TP , I X, PCS, PC NS , I DU M , NS , A VG ( 2 ) ,T AVGI 1000 ) , 
*NC ,TBAR,TSC,NGAME,NREPS ,NGRPS, I CUT ,OVRLAP , X3AR, SD 

COMPUTE THE AVERAGE NUMBER OF UMTS CORRECTLY 
CLASSIFING ThE CONTACT 

SLM=C. 0 
M = NOS I K ) 

CC 300 1=1 . M 

IF (RN2I IX l-PCNS) •‘.01,401,300 

401 S LM= SUM + 1 . C 
3 CC CONTINUE 

SL M = SUM /M 
L = 0 

IF THE AVERAGE IS GREATER THAN THE INPUT 
FRCBABIl ITY THE CLASSIFICATION IS CORRECT 

IF (SUN-PCNS) 4C3,4C3,4C2 

402 L= 1 

403 PFTURN 
END 



C 



C 

C 



SLGPCUT INE OVER!*) 



COMMON /SS/ XS( 10001 ,YS( 1000) ,VXS( 1000) ,VYS( 1000) , 

I ST ATI 1000 ) , TS I 1000 ) , VMAX,NUS( 1C00 ) ,NUM,MAX 
COMMON /HS/ XHI6) ,YH(6) , ICCDEI6) , RS , DELT , NFS , T L , TCS 
COMMON /UTIL/TP,IX,PCS,PCNS,IDUM,NS,AVG(2) ,TAVG( 1000) , 
te’T^SjIKfNGAFE.NREPS.NGRPS, I CUT , QVRL AP , X3 AR , SD 
IF ( ISTATI 1) ) 4CI,4C1 ,400 



IF THE SS WAS CETECTEC INCREMENT 
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TFt CC UT. TE K s f. STORE THfr TlMc 



-CC N S =N S+ 1 
N C = NC+1 

TA VG(NC ) =T P-TL 

IF THIS WAS THt- LAST KEPLICATICN CUMPUTE 
TFt AVFFAlt NUMBER OF DETECTIONS L RESET 
THt CCU NT IRS 0 R AN PC N NUMBER GENERATORS 

AG 1 IF ( N.G A t / N R F P S * N R b P S - N G A M F ) 50C,4CJ,50C 
» C 2 I N P = NG A M F / N U PS 

A V(} < I N [' ) = F I CAT! NS) /NRfPS 
C L M M V = R Nl ( C ) 

i:lmmy=rn?( c ) 

I X =- 1 X 
NS = 0 

IF THIS V.Ab THE LAST GAME CCMPIJTS 
STATISTICS 5 SIGNAL OUTPUT 

IF ( NC AMF-NRFPS^NGRPS ) 500,403,403 

403 CALL Mf MFI T ( AVG.NGRPS , XB Ak , SO ) 

IF (NC) 4C4, 404, 405 

4 C 4 T P AP =C* C 
T S C = 0 . 0 
CLIP M)c 

4 C c CALL MOMENT! r A VC , NC , TB AR , TSO) 

404 RETURN 1 
5 C C RETURN 

t N C 



SU’RPUTINE RHSUNEOO 

CCMMCN /SS/ XS(1000»,YS( 1000), VXS( 1C00) ,VYS< ICGO) , 

* I STAT( 1 000 ) , TS ( 1000 ) , VM AX , NUS ( 1 000 ) ,NtJM, MAX 
CCMMCN /HS/ XIU c ) , YH( 6) , ICCOE ( 6 ) , R S , DEL T , NH 5 , TL ,TCS 
CCMMCN /UT IL/TP, IX, PCS, PCNS, I DU M , N S , A VG < ? > , T A V G ( 1000) , 
*NC ,TBAR ,TSL ,NGAMh ,NRfPS, NGHPSt ICUT, CVRLAP ,XHAR, SB 

SET NUMBER CN THE SC t NE TC 0 4NC CANCEL 
Af v STORED EVENTS FCR THIS CENT ACT 

N C S ( K ) = G 

CALL RHM'LVF ( 3 , K ) 

CC 300 1=1 , NHS 

CALL bPKDIP FTP ALL UMTS CN THE SCENE 

IF (ICCCL(I)-K) 300,401,300 
4 C 1 CALL 6 K K T I P ( I ) 

3 C C CONTINUE 

r. c LFTE Thr FALSE CCNTACT 

I ST AT ( K ) =1 

RETURN 

ENC 



slhrcutine SEE 

CCMMCN /SS/ XS < 1000 I , YS ( 1000 ) ,VXS( 1000) , VYS( 1CCC) , 

* I STAT( 1 000 ) , TS ( 100C ) , VMAX ,NHS (1000 > ,NUM, MAX 
CCMMt.M /HS/ XH( 6 ) , YH( U ) , ICC06 ( 6 ) ,RS ,OELT ,NHS , TL, TCS 
CCMMCN /UT IL/TP , IX , PCS, PC NS, I DUN , NS , AVG ( 2 ) , T A VG ( 1 COC ) , 
*N( ,TBAR ,TSC,NGAME,NREPS,NGRPS, I CUT , CVRL AP , XG AR , SD 
RR ITE ( I PUT , 200 ) TP 

PRINT EACH HS PLSITICN L RANGE TO THE SS 
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H = SQP T ( ( XH (I > - X S ( 1 ) ) **2 + ( VH ( I )-YS( 1 » > > *2000.0 

VvR I T E ( I CUT , 202 ) I » XH ( I > , YH( I > » I CODE ( I) , P 
3C0 CCMINUE 

Ik R I TE ( I CUT ,201 ) 

FFINT EACH FALSE CONTACT POSITION £ CODE 

£f IT E ( I CUT * 202 ) I,XS ( I ) , YS ( I) » ISTAT ( I) 

3 C 1 CONTINUE 

2C0 FORMAT ( III 10X , 1 4HS T ATUS AT T I ME , F 8 . 3 , / / 20 X , 9HHS STATUS 
*// ) 

201 FORM AT ( 1H0 » 20X , 14HCGNTACT STATUS) 

2C2 FORMAT (1H0»I12»2F30.4» I12,F11»2 ) 

END 



SUBROUTINE SETUP(N) 

COMMON /SS/ XS ( 1000 ),YS( 1000 ) ,VXS( 1000) ,VYS( 1C00) , 

* I ST AT (1000) , TS( 1000), VM AX, NOS (1000) ,NUM,MAX 
COMMON / HS/ XH ( 6 ) , YH ( 6 ) , I C COE ( 6 ) ,RS,DELT,NHS,TL,TCS 
COMMON /UT IL/TP,IX,PCS,PCNS, I DU M , N S , A VG ( ? ) ,TAVG( 1000) , 
*NC ,TBAR ,TSC ,NGAMF,NREPS,NGRPS, I C UT , CVRL AP , X B AR , S D 
IF (N) 401 , 401 » 400 



IF N > 0 COMPUTE TIME INCREMENT FOR 
EACH CALL TD SEE £ STORE THE EVENT 



-,CC T1NC = ( TGS-TL) /n 
T 0 E = T L 

CC 300 1=1, N 

CALL SNfc( TOE, 6 , IDUM , IDUM , £999 ) 
I CE=TCF+T INC 
30C CONTINUE 
401 RETURN 
199 CALL C A L F U L ( 9 ) 

END 



SLERLUTINE TRACK(K) 

COMMON /SS/ XS( 1000) ,YS( 1C00 ) ,VXS( 100C) ,VYS( 1000) , 

* I ST AT ( 1000 ) ,TS ( 1000 ) , VMAX ,NOS ( 1000 ) ,NUM,MAX 
CCMMCN /HS/ XH(6> ,YH(6) , ICCD6(6) , RS , DFLT , NFS , TL , TOS 
COMMON /UTIL/TP,IX,PCS,PCNS, IDUM, NS, A VG (2) ,TAVG( 1000) , 
*NC ,TGAP , TSC ,NGAME, NRFPS, NGRPS, I OUT ,OVRLAP , X3AR , SD 
CALL UPDATE 

GET JOINT CLASSIFICATION 

CALL JOINT (K ,L ) 

IF (L) 401,401,400 

CONTACT HAS BEEN CCRRFCTLY CLASSIFIED 

4 C C CALL RESUMC(K) 

GOTO 503 

CONTACT HAS NOT BEEN CCRRFCTLY CLASSIFIED 
STCRF TFF NEXT TRACK EVENT £ CHECK THE NUMBER 
CN THE SCENE 

<,01 CALL SNE (TP + CELT ,4, I0UM,K, £999) 

IF ( NOS ( K) -MAX ) 403,502,502 

IF A HS IS AVAILABLE SEND HIM TC THE SCENE 

4 C 2 CC 300 1 = 1 , NHS 

IF ( ICODE ( I ) ) 5C0, 3C0, 300 
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JCC CCNTIMJF 

LF 301 1=1 , N h S 

I E ( ICOCf ( I ) ) 3C 1 , 5 CO » ? C 1 
? C 1 CONTINUE 
CC TO 302 

CC C/'LL OF MV t (2,1) 

Cl ICCCfc ( I ) = K 

N C S ( K ) = NCS { N ) + 1 

LOCATE TFE POSITION CF ALL HS CN THE SCENE 

C 2 CC 302 1=1 , NHS 

IF ( ICOFCI ! )-K) 3 C 2 , * ♦ C 2 ,302 

02 XI- < I ) = X S ( K ) 

YE ( I ) = YS (K ) 

C 2 CONTINUE 

03 PC TURN 

99 CALL CAl FUI. ( A ) 

:n c 



SLEMLTINE UPDATE 
C 

CCM'CN /SS/ XS( 1000), YS( 1000), VXS( 1C00) ,VYS( 1C00) , 

* I S T A T { 1 000 ) , TS ( 100 0 I , V P A X , NO S ( LOGO ( , NUM , MAX 
CCPMN /HS/ XH ( 6 ) , YH( fc ) , ICCOC ( fc ) , RS ,l)tL T ,NH.S , Tl , TCJS 
C( P P, C N /UT 1L/TP, IX, PCS, PCNS, I OOP , NS , AVG ( 2 I ,TAVG( 10001 , 
MC , TCAR , TSC , NGAPC , NMPS , NGMS, I CUT , CVM AP ,X3 Ap f SC 
CL 3CC I - 1 , N U M 
C 

C IE TF-C CONTACT IS ACTIVE UPCATE ITS POSITION 

C 

IF ( 1ST AT ( I ) — 1 I ^OU, 300,400 
MC XS(I)=XS(I ) + V X S ( I ) M TP-TS( I ) ) / 6 C . 0 
Y S ( I ) = YS( I ) +VYS( I ) M T P - T S < I ) ) / 60. 0 
T $ ( 1 ) = T P 
3CC CLNTINU C 
FCTUPN 
ENC 
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